DPNU-96-40 
August 1996 



On Existence of Nontrivial Fixed Points 
in Large N Gauge Theory in More than Four Dimensions 



Jun NishimuraJ] 

Department of Physics, Nagoya University, 
Chikusa-ku, Nagoya 4-64-01, Japan 



Abstract 

Inspired by a possible relation between large N gauge theory and string theory, we 
search for nontrivial fixed points in large N gauge theory in more than four dimensions. 
We study large N gauge theory through Monte Carlo simulation of the twisted Eguchi- 
Kawai model in six dimensions as well as in four dimensions. The phase diagram of the 
system with the two coupling constants which correspond to the standard plaquette 
action and the adjoint term has been explored. 
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1 Introduction 



String theory has been considered as a natural candidate for the unified theory including 
gravity. Intensive study revealed, however, that nonperturbative understanding of its dy- 
namics is essential in making any predictions as to our present world. The double scaling 
limit of the large N matrix model is a rare example in which nonperturbative study 
of string theory has been successful. Although the development concerning matrix models 
so far is restricted to the space-time dimension less than 1, the above success shows that 
nonperturbative study of string theory in physical dimensions may be possible through the 
study of the double scaling limit of some kind of large N field theory. 

As such a field theory, we consider large N gauge theory. It has been generally believed 
that large N gauge theory is related to some kind of string theory, since 't Hooft showed 
within perturbation theory that the 1/N expansion of U(iV) gauge theory gives an expansion 
in terms of the topology of the Feynman diagrams 0. If, in the strong coupling regime, 
dominant contributions come from such diagrams as those with very fine internal lines like 
fishnets, the theory can be viewed as a string theory by identifying the fishnet with a con- 
tinuous worldsheet embedded in the space-time where the gauge field lives. Also in lattice 
gauge theory, the expectation value of a Wilson loop can be expressed, in the strong cou- 
pling expansion, as a sum over random surfaces spanning the loop. Although no rigorous 
connection between large N gauge theory and string theory has not yet been established, 
there is a good reason to consider it seriously. 

In four dimensions, since nonabelian gauge theory is asymptotically free, a possible rela- 
tion to string theory can arise only as an effective theory in the infrared region, where the 
coupling constant becomes sufficiently large. This naturally leads us to consider the theory 
in more than four dimensions. 

The e-expansion of (4 + e) -dimensional nonabelian gauge theory shows that there exists 
an ultraviolet fixed point between the confining phase and the deconfining phase when e is 
positive. If this fixed point survives even in physical dimensions such as 5 or 6, there is a 
possibility of constructing a string theory by approaching the fixed point as one take the 
large N limit. Considering that the critical dimension of bosonic string theory is 26, one 
may speculate that the fixed point might well exist in more than four dimensions. 

Since gauge theory is perturbatively unrenormalizable in more than four dimensions, 
the existence of such a fixed point is quite nontrivial. An example is three-dimensional 



nonlinear sigma model, which is perturbatively unrenormalizable and yet can be defined 
on the lattice by taking the continuum limit at the second order phase transition point. 
Although this example is rather trivial in the sense that the continuum theory belongs to 
the same universality class as three-dimensional 4 theory which is superrenormalizable, it 
shows that constructive definition of perturbatively unrenormalizable theory is possible at 
least in principle. Constructing quantum gravity within ordinary field theory belongs to the 
same kind of problem and we expect that our study may also be suggestive to that direction. 

Since we need a nonperturbative approach, we consider lattice gauge theory. Monte Carlo 
simulation of large N gauge theory in a straightforward manner is, however, unfeasible due 
to the huge number of dynamical degrees of freedom, which is as large as iV 2 x D x L D , where 
L is the size of the lattice, and D is the dimension of the space-time. However, as Eguchi 
and Kawai showed in the early 80's 0, the factorization valid in large N gauge theory leads 
to the fact that the above system is equivalent to a one-site model {L=l) with N 2 x D 
dynamical degrees of freedom. This tremendous reduction of dynamical degrees of freedom 
makes the study of large N gauge theory accessible by numerical simulation. 

This paper is organized as follows. In section 2, we explain the reduction of Wilson's 
lattice gauge theory to the one-site model in the large ./V limit. In section 3, we study the 
four- dimensional case. In addition to the standard plaquette action, we put the adjoint term 
in the action and study the system with the two coupling constants. In section 5, we show 
our results for the six-dimensional case. Section 6 is devoted to conclusion and discussion. 

2 Reduction to one-site model 

In Wilson's lattice gauge theory, SU(iV) gauge theory can be defined with the action 

S = -N(3J2 E tr(U n> ,U n+ ^Ui + ^Ul,v)- (!) 

?7 n . M 's are SU(iV) matrices and are called link variables. A gauge invariant observable can 
be given by the Wilson loop 

W(C) = ^tr(U nja U n+&tP U n+&+ ^ ■ ■ ■ U^), (2) 

which is the trace of the ordered product of link variables along the loop C. 

The one-site model which is equivalent to the above Wilson's theory can be defined with 



the action 

S = -N/3 Y, Z,MU,UMlUi), (3) 

where Z^ v is an element of the center group Zn and is called 'twist'. A successful choice of 
the twist is given by 

Z^y = exp(27ri/L) for ji < u, (4) 

%vn = Z* u , (5) 

where L is an integer defined through 

N = L D ' 2 (6) 

with the space-time dimension D, which is taken to be even []. The twist has been introduced 
by the authors of Ref . |4[] in order to cure a problem in the original Eguchi-Kawai model || , 
where Z^s were taken to be unity. The model is thus called twisted Eguchi-Kawai model. 
For a general review on twisted Eguchi-Kawai models, we refer the reader to Ref. || . 

The observable corresponding to the Wilson loop (|2]) can be defined in the twisted Eguchi- 
Kawai model as 

W ( C ) = 4(11 Z^»MUaUpU r ■ .U u ), (7) 

where Np^ is the number of plaquettes in the {n,v) direction on the minimal surface spanning 
the loop C. 

Under the assumption of the factorization 

(W(Ci)W(C 2 ) • • • W(C k )) = {W{C l )){W{C 2 )) ■ ■ ■ (W(C k )) + O (J^j , (8) 

which is valid generically in large gauge theory, one can show that (w(C)) calculated 
in the twisted Eguchi-Kawai model is equal to (W(C)) calculated in the original Wilson's 
theory on the LP lattice, where L is related to iV through eq. (^). 

In this sense, the twisted Eguchi-Kawai model is equivalent to the original Wilson's theory 
on the infinite lattice, in the large iV limit. As is indicated in the above statement, the finite 
iV effects in the twisted Eguchi-Kawai model appear in two ways. One is the violation of 
the factorization which is the assumption of the equivalence, and the other is the finite L 
effects in the corresponding Wilson's theory to which the one-site model is equivalent. 

1 For odd D, one can apply the reduction scheme only to the (D — 1) dimensions leaving one dimension 
unreduced. We prefer to restrict ourselves to even dimensions, i.e. D — 4 and D = 6 in order to avoid any 
technical complications. 



Let us now explain how to perform Monte Carlo simulation of the twisted Eguchi-Kawai 
model. In contrast to the ordinary Wilson's theory, the action (^) is not linear in terms of 
each link variable. In order to make the heat bath algorithm || applicable to the model, 
we use the technique proposed by Ref. @. The idea is to introduce an auxiliary field 
Qfiu{^ < A* < v ^ D), which is a general complex N x N matrix, with the following action. 



s = NpY.trQtrQr 

\i<V 

-N/3 trQl^U^ + t vt p v V^ 
-N(3 £ trQ^Mul + tlpUpl), 



fl<U 

where t^ u is the square root of Z^ v . Since this action is linear in terms of U^, we can use 
the heat bath algorithm for the update of U^. We update by successively multiplying 
it by matrices each belonging to the N(N — l)/2 SU(2) subgroups of the SU(iV) ||. After 
updating all the U^s in this way, we perform the update of Q^, which can be done with 
little cost by generating Gaussian variables. This defines the 'one sweep' of our system. For 
further technical details of the algorithm, we refer the reader to Ref. 0. 

3 Results for four-dimensional case 

We first study the model in four dimensions. We take N = 16 which corresponds to L = 4 
in the equivalent Wilson's theory. In fig. 1, we plot the mean plaquette defined by 

P = (^ d{ d- 1) Y,/»MU»U v UtU^ . (9) 

Each point is an average over 1000 sweeps. We also plot the leading terms of the strong- 
and weak-coupling expansion for the Wilson's theory. 

j (for small 0) 

\^r P (for large 0). ^ 10 > 

Our data are in good agreement with those of Ref. J|]. We observe a clear indication for a 
first order phase transition at (3 ~ 0.35. We should note that a first order phase transition 
is already present at finite N, for N > 4, when one uses the standard plaquette action 0. 
This first order phase transition is considered as a lattice artifact of the plaquette action, 



and can be removed, for example, by using the modified action 

S = -NPj2trU(P) - ^£tr A [/(P), (11) 
p z p 

where U(P) is the ordered product of link variables around the given oriented plaquette 
P. The second term is the trace of U(P) taken in the adjoint representation, which can be 
rewritten as 

tr A f/(P) = |trf/(P)| 2 -l. (12) 



Indeed, in Ref. |Tl| , Monte Carlo simulation of the above system in the ordinary Wilson's 
theory has been done for 4 < iV < 8 and it is shown that the line of first order phase 
transition terminates at a point with a sufficiently large negative (3 a- 

We, therefore, study large N gauge theory with this modified action. The adjoint term 
can be incorporated into the twisted Eguchi-Kawai model with the action 

S = -N/3 £ Z^MU,U u UlUl) - y E \Z^iU,U v UlUt\ 2 . (13) 

The second term, which corresponds to the adjoint term, can be dealt with by Metropolis 
algorithm. In order to increase the acceptance rate, we rewrite the action as 

S = -N((3 + (3 A C) £ Z^MU,U v UlUl) - y E I Z„ v txUJJ u UlUl - NC\ 2 , (14) 

taking the arbitrary real constant C to be approximately equal to the mean plaquette defined 
byeq. (|). 

In fig. 2, we show the mean plaquette as a function of (3 for (3 A = —0.3, —0.7, —1.0. Each 
point is an average over 1000 sweeps. We replot the data for f3 A = 0. As Pa goes to a 
larger negative value, the gap at the transition point decreases, and finally disappears. This 
implies that the first order phase transition can be removed by introducing a sufficiently 
large negative fix even in the large N limit. 

It is known that, in the large N limit, the theory with ((3, (3a) is equivalent to the theory 
with (/?', (3' A ), if 

P> = P + (Pa-P'a)P((3,{3a), (15) 
where P(/3,/3a) is the mean plaquette at the coupling constants ((3, (3a) |~2], [13 . 



2 There may be various types of modifications that can be used instead of ([tl|). In Ref. JHJ, it is reported 
that the first order phase transition is absent when one uses the renormalization-group improved action. 



Using this equivalence theorem, we can predict the results for /5a = with the input 
of the data for /3a = —0.3, —0.7, —1.0. In fig. 3 we show the prediction from the data 
for Pa = —0.3,-0.7,-1.0 together with the data for /3a = 0. Note that there are three 
predicted values for each /3 in the critical region. The middle point is unstable and can never 
be seen in the simulation. The other two points are stable or metastable. The metastable 
vacuum is stabilized as N increases due to the small tunneling probability ~ exp(— const. N 2 ). 
Each of the three points corresponds to a solution to the Schwinger- Dyson equation for the 
Wilson loops fll3"| . The 'equivalence' of the systems with different /3 A 's should, therefore, be 



understood in the sense that the two systems have a common solution to the Schwinger- 
Dyson equation. One can see that our data clearly satisfies the equivalence theorem. 

4 Results for six-dimensional case 

Let us turn to the six-dimensional case. In fig. 4, we show the mean plaquette as a function 
of /3 for N = 64 (L = 4) with the standard plaquette action (/3a = 0). Each point is an 
average over 300 sweeps. We see a strong hysteresis indicating a first order phase transition. 
Note that the thermal cycle is substantially larger than that in the four-dimensional case. 
Fig. 5 shows the result for N = 27 (L = 3) with /3a = —10.0. Each point is an average over 
300 sweeps. Here we see a striking difference from the four-dimensional case. We see that 
there are two different /3's that give the same value for the mean plaquette. It is obvious 



from the relation (15) that the phase transition will never become continuous, however large 
a negative value of /3a we may take. 

Using the equivalence theorem in the previous section, we can examine the asymptotic 
behavior of the system in the large negative /3a region. From the relation fllBD , we have 

P(/3,/3 A ) = P(/3+ (/3 A -/3 a )P(/3,/3 a ),/3 a ). (16) 

We fix j3 and /3a and consider the asymptotic behavior of the system with /?' = /? + (/?a — 
/3 A )P(/3, /3a) and f3' A for j3' A — > — oo. Differentiating the above equation with respect to /3, 
we obtain 



dP 
~df3 



' (17) 



Therefore, denoting the horizontal distance of the two nearly parallel data lines in fig. 5 as 



Aft, we can write the gap in the mean plaquette as 

dP 



AP = Aft • 



dft 



Aft_ 



(18) 



Also the fluctuation of the mean plaquette P can be written as 
5P = 



(p 2 ) - (py< 



1 dP 
Jp ~dft 



Wa\ N 



(19) 



Thus the fluctuation SP as well as the gap AP vanishes in the ft' A — > — oo limit. Note 
also that the gap AP vanishes faster than the fluctuation SP, so that the discontinuity, 
which is a clear signal of the first order phase transition, becomes difficult to observe in 



the actual simulation. These results naturally explain the observations of Ref . |TJ[ , where 
five-dimensional pure SU(2) lattice gauge theory has been studied. 

In the large N limit, by virtue of the equivalence theorem, we can conclude that the line 
of first order phase transition continues to (ft, /?a)=(oo, — oo). Considering that the data of 



Ref. [JL4J] are qualitatively the same as ours, we may naturally expect that the same is true 
for 2 < iV < oo. 

One may think of enlarging the coupling-constant space where one searches for a second 
order phase transition. For this purpose, we consider those additional terms in the action 



which can be obtained from (jTT|) by replacing plaquettes with three-dimensional loops which 
consist of six links. Such terms can be incorporated into the twisted Eguchi-Kawai model 
by the following additional terms. 

~7a£ E \Z» u Z uX Z^tr(U,U v UxUlUtu{) - Na\ 2 

-7a E \z; u z: x zi^(u,uiu x ulu u u{)-Na\ 2 . 

jjL<V<\ 

Corresponding to the mean plaquette, we define the mean three-dimensional loop by 



Q 



l 

N 



D(D-1)(D-2) ^ X ^^{ Z ^ Z lM U ^ U ^ U t U l)) 



6 



£>(£>- !)(£>- 2) 



E Re{z; v z: x zi^(u,uiu x ulu u ui)} 



li<v<\ 



A similar argument as with ft, ft a and P applies to 7,7a and Q, where 7 = 7a«- 

In fig. 6, we show the data for (P,Q) in the ordered phase and in the disordered phase 
for four sets of (ft,j) with (3^ = —1.0 and 7a = —2.0. Each point is an average over 50 



sweeps. Note that the quadrilaterals spanned by the data points in the two phases have an 
overlap. This means that the same (P,Q) can be obtained for different (/3,7)'s with fixed 
ft a and 7a- We can thus conclude that the phase transition will not become continuous even 
in the enlarged coupling- const ant space. 

5 Summary and Discussion 

In this paper we study large N gauge theory in six dimensions as well as in four dimensions. 
The phase diagram in the large negative ft a region has been clarified by making use of a 
simple property of the system with the two coupling constants ft and Pa due to the large N. 

In four dimensions the line of first order phase transition terminates at a sufficiently 
large negative ft a and the two phases are actually connected analytically. At the end point 
of the line, the first order phase transition becomes second order. What is the continuum 
theory that can be defined at this end point? Since the phase transition is not a deconfining 
one, the string tension does not scale to zero at the end point. This means that minimum 
surfaces dominate the summation over surfaces spanning the Wilson loops. Considering, for 
example, the two-point function of plaquettes, dominant contributions come from tube-like 
surfaces connecting the two plaquettes. This means that the continuum theory looks more 
like particle theory rather than string theory. This is consistent with the claim for SU(3) 
gauge theory that the continuum theory defined at this end point is (trivial) <j) A theory [ |15|] . 

In six dimensions, the line of first order phase transition continues to infinity. Moreover, 
we see that this is the case even if we enlarge the coupling-constant space to include three- 
dimensional loops in addition to plaquettes. It is natural to expect that this pattern is 
repeated for other types of loops. Thus we consider that there is no nontrivial fixed point 
in large N gauge theory in six dimensions. 

In the context of string theory, it is tempting to interpret our conclusion as a result 
of tachyon instability of bosonic strings, which may be cured by introducing space-time 
super symmetry. It is, therefore, probable that we may find a nontrivial fixed point in six 
dimensional super symmetric large N gauge theory, though it will require much more effort 
since we have to deal with dynamical fermions. 
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Figure captions 



Fig. 1 The mean plaquette P as a function of the coupling constant (3 in four-dimensional 
SU(16) gauge theory with the standard plaquette action ((3a = 0). The dashed line represents 
the strong coupling expansion P — (3, whereas the dash-dotted line represents the weak- 
coupling expansion P = 1/(8(3). 

Fig. 2 The mean plaquette P as a function of the coupling constant (3 in four- dimensional 
SU(16) gauge theory with (3a = —0.3 (triangles),— 0.7 (squares),— 1.0 (diamonds). We replot 
the data for (3a = (crosses) shown in fig. I. 

Fig. 3 The predicted values for the mean plaquette P as a function of the coupling constant 
(3 for (3a = with the input of the data for (3a = —0.3 (triangles), —0.7 (squares), — 1.0 
(diamonds) shown in fig. 2. We also plot the original data for (3a = (crosses). 

Fig. 4 The mean plaquette P as a function of the coupling constant (3 in six-dimensional 
SU(64) gauge theory with the standard plaquette action ((3a = 0). The dashed line represents 
the strong coupling expansion P — (3, whereas the dash-dotted line represents the weak- 
coupling expansion P = 1/(12(3). 

Fig. 5 The mean plaquette P as a function of the coupling constant (3 in six-dimensional 
SU(27) gauge theory with (3a = —10.0. The triangles are the data in the disordered phase, 
whereas the squares are the data in the ordered phase. 

Fig. 6 The mean plaquette P and the mean three-dimensional loop Q in six-dimensional 
SU(27) gauge theory are plotted for four sets of ((3,^) with (3a = —1.0 and 7a = —2.0. 
The triangles are the data in the disordered phase, whereas the squares are the data in the 
ordered phase. 
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